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ABSTRACT 

Context. Multiwavelength observations have revealed the highly unusual properties of the y-ray source PMN11603-4904, which are 
difficult to reconcile with any other well established y-ray source class. The object is either a very atypical blazar or compact jet 
source seen at a larger angle to the line of sight. 

Aims. In order to determine the physical origin of the high-energy emission processes in PMN 11603-4904, we study the X-ray 
spectrum in detail. 

Methods. We performed quasi-simultaneous X-ray observations with XMM-Newton and Suzaku in 2013 September, resulting in the 
first high signal-to-noise X-ray spectrum of this source. 

Results. The 2-10keV X-ray spectrum can be well described by an absorbed power law with an emission line at 5.44 ± 0.05 keV 
(observed frame). Interpreting this feature as a Ka line from neutral iron, we determine the redshift of PMN 11603-4904 to be 
z = 0.18 ± 0.01, corresponding to a luminosity distance of 872 ± 54Mpc. 

Conclusions. The detection of a redshifted X-ray emission line further challenges the original BE Lac classification of 
PMN 11603-4904. This result suggests that the source is observed at a larger angle to the line of sight than expected for blazars, 
and thus the source would add to the elusive class of y-ray loud misaligned-jet objects, possibly a y-ray bright young radio galaxy. 

Key words. Galaxies: active - Galaxies: jets - Galaxies: individual: PMN 11603-4904 - X-rays: galaxies - Gamma rays: galaxies 


1. Introduction 

Active Galactic Nuclei (AGN) are among the brightest extra- 
galactic sources at X-ray energies. According to the unification 
model ( |Antonucci|1993[[Urry & Padovani|1995| l, the orientation 
of AGN with respect to the line of sight, different accretion rates, 
and the presence of a powerful jet lead to different spectral ap¬ 
pearances. Jet-dominated sources at the smallest jet inclination 
angles are classified as blazars. Sources which are oriented at 
larger angles to the line of sight are typically less dominated by 
beamed jet emission. The origin of the X-ray emission of radio- 
loud AGN is still an open question (e.g., |Fukazawa et al.|2014| l. 
It can arise from Comptonization of seed photons from the ac¬ 
cretion disc in the hot corona, synchrotron or inverse Compton 
emission, or a combination of these. The most prominent emis¬ 
sion line in AGN X-ray spectra is the Fe Ka line at ~6.4keV 
(e.g., |Nandra|2006] Shu et al.|2010 Ricci et al.|2014| l. This line 
is usually explained by fluorescence of neutral or mildly ionized 
iron produced by the primary, hard X-ray photons irradiating the 


accretion disk or the molecular torus. Measurements of iron line 
emission probe therefore the inner regions of AGN. The detec¬ 
tion of an iron line is a clear indication of reprocessed radiation 
from matter in the vicinity of the primary X-ray source. Iron 
lines are common features in spectra of radio-quiet and radio- 
loud sources at larger jet angles to the line of sight ( |Bianchi 
et al.|2004l. In contrast, typical blazar X-ray spectra are feature¬ 


less continua, clearly dominated by the beamed jet emission that 
outshines any putative underlying emission components (e.g., 
[Rivers et al.|20I3]|Grandi et al.|200^ . 


In jMiiller et al.| ( |20I4| hereafter Muel4), we discussed the 
multiwavelength properties of PMN J1603-4904, which is asso¬ 
ciated with a bright y-ray source with a remarkably hard spec¬ 
trum (2FGL 11603.8^904; [NoTan et al.||2m2l l. These proper¬ 
ties are difficult to reconcile with its original classification as a 
BLLac type blazar ( Nolan et al.|2012||Ackermann et al.|2013| l: 
High resolution TANAMI ( Ojha et ail.||2010| l radio observa- 
tions using Very Long Baseline Interferometry (VLBI) revealed 
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a symmetric structure on milliarcsecond (mas) scales with the 
brightest, most compact component at the center. The broadband 
spectral energy distribution (SED) showed an excess in the in¬ 
frared, which could be modeled with a blackbody spectrum with 
T ~ 1600 K, consistent with the emission of circumnuclear dust 
heated by the disk photons. Such thermally reprocessed radia¬ 
tion is not expected in typical blazar SEDs. Optical measure¬ 
ments constrained the redshift to z < 4.24 ( |Shaw et al.||2013| l. 
Based on the results presented by Muel4, however, no conclu¬ 
sive classification of PMN J1603-4904 could be made. While 
Swift/XRT observations showed a faint X-ray counterpart, low 
photon statistics did not allow us to sufficiently constrain its X- 
ray spectrum. 

In this letter, we report on the first XMM-Newton and 
Suzaku observations of PMN J1603-4904 conducted within the 
scope of the TANAMI multiwavelength program (Sect. [^. In 
Sect. we present the X-ray spectral characteristics and discuss 
its implications in Sect. 


2. Observations and data reduction 


In 2013 September, we performed quasi-simultaneous X-ray 
observations of PMN J1603-4904 with XMM-Newton (|Jansen| 
[eTaLjlWTI obs-id 072470010 1, performed 2013-09-17) and 


Suzaku (Koyama et al. 2007 obs-id 708035010, performed 


2013-09-13) in order to constrain the X-ray spectrum better than 
previous observations with Swift/XRT ( [Burrows et al. 2005 1 
Muel4). The XMM-Newton data processing, analysis, and ex¬ 
traction of the spectrum were performed using the standard tasks 
of the XMM System Analysis Software (SAS 13.5.0). We de¬ 
tected a single, unresolved. X-ray source at the radio position 
of PM N J1603-4904. Only data from EPIC-PN ( jStriider et ah 
exposure time 39.0ks) and EPIC-MOS 1 (Turner et al. 
exposure time 33.4ks) cameras were used, as the usable 


2001 


2001 


exposure of EPIC-MOS 2 was only 4.7 ks. The PN spectrum was 
extracted using a circle with 35" radius centered on the coordi¬ 
nates of the radio source. Eor the background spectrum a circular 
source-free region of 50" radius was chosen. Eor the MOS 1 data 
extraction, background and source regions with radii of 100" 
were used. We chose larger extraction radii for the background 
to minimize the measurement uncertainty. We used the HEA- 
SOET v6.15 package for the Suzaku data analysis. Extraction 
and background regions of 94" were used for all XIS detec¬ 
tors, which after hltering had an exposure time of 48.5 ks. The 
spectral analysis was performed using the Interactive Spectral 
Interpretation System (ISIS, Version 1.6.2-27, [Houck & Deni-j 
|cola|2000) l. An independent analysis of the XMM-Newton and 
Suzaku data revealed no signihcant flux variability during the 
observed period. In order to improve the statistics and ensure the 
validity of the ;\f^-statistics, all 1-10 keV spectra were rebinned 
to a minimum signal-to-noise ratio of 5 and modeled simulta¬ 
neously considering the background. Throughout this paper we 
use the standard cosmol ogical model with Qm = 0.3, A = 0.7, 
Ho = 70kms“*Mpc“' ( Beringer et al.|20T2 i. All uncertainties 
quoted in the following are given at the 90% conhdence. 


3. Results 

3.1. X-Ray continuum 

The quasi-simultaneous data from XMM/PN, MOS 1, and 
Suzaku/XIS were simultaneously modeled with an absorbed 
power-law component, a Gaussian emission line, and cross¬ 
calibration constants to account for the relative flux calibrations 
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Fig. 1. Simultaneous fit to the XMM-Newton and Suzaku data (see Ta¬ 
ble [TJ. The data are best fitted by an absorbed power-law component 
with an emission line at ~ 5.44 ± 0.05 keV. (a): Count spectrum for all 
detectors fitted with an absorbed power law and a Gaussian emission 
line (models shown in gray). Ratio of data to model for a fit of an ab¬ 
sorbed power law (b) and including a Gaussian emission line (c). (d): 
Unfolded, combined spectrum of all data sets with the best-fit model 
(red). The shaded region highlights the position of the emission line. 


of the instruments. We use the tbnew model ( [Wilms et al.|2012| l 
to account for neutral Galactic a bsorption (Galactic H I value is 

Kalberla et al.||2005| , and the 


VI. 


cross sections and abundances of[Verner et al.|([1996[l and[Wilms[ 


et al. (2000 1 , respectively. Eurther residuals reveal the presence 


of source intrinsic absorption, which we model with an addi¬ 
tional, redshifted absorption component (tbnew_z, see below). 
We detect a strong emission line that can be modeled by a Gaus¬ 
sian component at 5.44 + 0.05 keV with an equivalent width of 
EW = 200 + 90eV (Eig.[T}i). Using the centroid energy of the 
Gaussian emission line with respect to the neutral Ee Ka line 
line rest frame energy, we can constrain the redshift of the sys¬ 
tem to z = 0.18 + 0.01 (see Sect. [3.2[ for details). We find a 
photon index of Tx = 2.07/° j 2 with an intrinsic absorption of 
Ah = 2.05/q[2 X lO^^cm^^. The best fit model parameters can 
be found in Table[T] Adding a Gaussian line to describe the emis¬ 
sion line improves the fit fmmx^ - 198.1 (165 degrees of free¬ 
dom [dof]) to - 183.3 (162 dof). To test for the significance 
of the line feature at 5.44 keV, we performed a Monte Carlo sim¬ 
ulation based on the Null-Hypothesis that the intrinsic model for 
the measured counts is an absorbed power-law lacking a line. 
This hypothesis is rejected with a p-value of p < 1.1 x 10 
Eigure [T] shows the spectra of all single datasets with the best 
fit model and the combined datasets regridded to a joint energy 
grid. 
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Table 1. Best fit parameters using an absorbed power-law model with a 
Gaussian emission line. 


Parameter 
A^h (intrinsic) 

■P' 2 -iokeV (de-absorbed) 

r 

Crpe Kff 

EW 

constpn 

constMosi 

constxiso 

constxisi 

constxis3 _ 

TTd^f 


Value 

2.05!“J^x 10^^ 

(4.39 + 0.17) X 10“'^ erg s“* 

9 A7+O.O4 
- 0.12 

5.44 +0.05 keV 
< 0.06 eV 
200 + 90eV 
1.0 

-0.08 

1.57 + 0.11 
1.46 + 0.12 
1.63 + 0.12 
183.0/162 


3.2. Iron emission line 


Due to its high fluorescent yield, the Fe Kcr transition produces 
the most prominent line feature in the 2-lOkeV energy range. 
Such lines are seen in the X-ray spectra of the majority of 
Seyfert and radio galaxies ( |de La Calle Perez et al. 2010 1 . In 
most sources where redshifts are known (e.g., Markowitz et al.| 
2007 1 Shu et al. 2010[ Ricci et al. 2014| l, the line energy and 
equivalent width are consistent with an origin by fluorescence 
with rest-frame energy of ~6.4 keV in a neutral to mildly ionized 
medium that is irradiated by hard X-rays. Even in sources show¬ 
ing relativistically broadened Fe Ka lines from the innermost 
regions of the accretion disk, narrow line components consistent 


with a rest-frame energy of 6.4 keV are typically seen (Risaliti 
|et al.|[20T3) |Marinucci et al.|[2014a|b| l. We therefore consider it 
most probable that the emission line seen in PMNJ1603-4904 
at 5.44 + 0.05 keV is also emitted intrinsically at 6. 4ke\Q such 
that the corresponding redshift of the source can be constrained 
to z = 0.18 + 0.01. 

Adopting this redshift, the luminosity distance is 872 + 
54Mpc. In the rest-frame, the intrinsic X-ray luminos¬ 
ity, neglecting Galactic and source intrinsic absorption, is 
E2-iokeV = (4.03 + 0.20) X 10"^^ergs“' and the y-ray lu¬ 
minosity is Li-iooGeV = (8.1 + 1.1) X 1 0"^^ergs~* (for 

S 1-100 Gev = (1.44 + 0.07) X 10 


erg s ' cm 


Nolan et al. 


2012)1. At 8.4 GHz, we derive a rest-frame luminosity density 

- ' i__ 1,- . . ~ 


of L8.4GHZ = (7.3 + 1.5) X 10^ ergs“ Hz“ (based on the flux 
from Muel4). The comparison of ATCA observations at arcsec 
resolution and VLBI observations with milliarcsecond resolu¬ 
tion shows that ~80% of the radio emission is emitted on scales 
of ~15mas (Muel4), i.e., ~46pc at z = 0.18. The remaining 
~20% is emitted on larger scales which are constrained to < 1" 
(Muel4), i.e., ~3kpc atz = 0.18. 


4. Discussion 


The X-ray properties derived from the XMM-Newton and 
Suzaku spectra allow us to further constrain the classification of 
this unusual y-ray source. Here we compare the X-ray spectrum 
of PMNJ1603-4904 to the typical spectra of different radio- 
loud AGN (see Table]^. 

Compared to other source types, blazars are the most lumi¬ 
nous AGN at the X-ray energies (e.g., Chang|2()T0 Rivers et al. 


* In the very unlikely event that the line is instead due to emission from 
a strongly ionized plasma, for hydrogen-like iron, the Lya resonance at 
6.966 keV would yield a redshift of z = 0.28 ± 0.01. 


2013 1 . With an X-ray luminosity of T 2 -iokeV ~ 4 x 10"^^ erg s“', 
PMN J1603-4904 is only comparable to the fainter blazars of 
the BLLac type. Its y-ray luminosity of Ly ~ 8 x lO'^^ergs^' 
is comparable to low or interm ediate peaketp] BE Lac objects 
and flat spectrum radio quasars ( |Ackermann et al.|201 l| l. Since 
blazar X-ray spectra are dominated by beamed jet emission, they 
are generally featureles^and can in most cases be modeled with 
a simple or broken power law ( |Chang|2010[ |Ushio et al.|2010) 
Rivers et al.|2013)l. The strong iron line in the X-ray spectrum of 


PMN J1603-4904 therefore suggests a non-blazar nature. The 
intrinsic X-ray absorption (Ah ~ 2x10^^ cm^) is furthermore un 
usually high compared to other blazars (Kubo et al. 1998) Chang 


2010) Rivers et al.|2013 1. We therefore conclude that its original 


classification as a BLLac ( |Nolan et al.|2012)|Shaw et al.|2013l l 
is very unlikely. 

The multiwavelength properties of PMN J1603-4904 
opened the door to an alternative classification as a jet source 
seen at a larger inclination angle, i.e., a radio galaxy (Muel4). 
Overall, the X-ray properties of PMN J1603-4904 are com¬ 
parable to those of radio galaxies (Table |^. In about 50% 
of all broad-line radio galaxies, Fe Ka lines with equivalent 


widths of typically a few 100 eV are detected (e.g., Sambruna 
|et al.||T999) . PMN J1603-4904 has a comparable equivalent 
width of EW = 200 + 90 eV. We caution, however, that only 


few non-blazar sources are detected at y-ray energies (Abdo 


|et al.|[2OT0t |Nolan et al.||2012| |Katsuta et al.||20l3] ). The y-ray 
luminosities of such “misaligned” jets show a broad range from 
Ly ~ 10^*-10"^®ergs“\ with FR I galaxies being typically less 
y-ray luminous than FR II sources. The y-ray luminosity of 
PMN J1603-4904 is only comparable to the most powerful 
misaligned sources ( Abdo et al.|2010| l. 

Among the sources seen at a large inclination angle. Com¬ 
pact Steep Spectrum (CSS) and Gigahertz Peaked Spectrum 
(GPS) sources can be interpreted as the younger versions of 
evolved radio galaxies, often showing compact symmetric ob¬ 
ject (CSO) radio morphologies at parsec scales (Readhead 
|et al.| |1996a|b[ |0’Dea |1998| l. Owing to its milliarcsecond- 


scale properties, we discussed the possible classification of 
PMNJ1603-4904 as a CSO in Muel4. The 2-10keV spec¬ 
tra of GPS and CSO sources can generally be well modeled 


by absorbed power laws (e.g.. 

Vink et ai.||2006[ jGuainazzi) 

et al.||2006[ Siemiginowska et al. 

|2008 

Tengstrand et al. 2009) 

Kunert-Bajraszewska et al.||2014 

Fe 

Ca emission is only de- 

tected in a few sources (e.g.. 

Guainazzi et al.||2006| |Siemigi-| 

nowska||2009| jTengstrand et a 

.12009 1 . 

The X-ray properties of 


PMN J1603-4904 are a good match to the findings in these sam¬ 
ple studies of GPS and CSS sources. These sources generally 
have a radio luminosity density of L > lO^^-lO^® erg s ' Hz ' 


at 5 GHz ( |0’Dea & Baum |1997| for combined complete sam¬ 
ples). Thus PMN J1603-4904 is compatible with the less pow¬ 
erful GPS and CSS sources. Its linear size is comparable to the 
canonical size limit for CSOs (< 1 kpc) or medium-size symmet 
ric objects (< 15 kpc; Fanti et al.|1995[ Readhead||1995[ 


Read- 


head et al. 1996b I. Gamma-ray emission was predictec 
theoretical modeling of GPS sources (e.g.,|Stawarz et alT 


from 

2008) 


^ The synchrotron peak frequency derived from the parametrization 
of the broadband spectral energy distribution is Vsy„c - 2.2 x lO'^ Hz 
(Muel4), which would be more in line with a low-peaked obje ct (de 
fined by < lO'"* Hz; 


peak 

As a rare exception. 


Fossati et al. 


1998 


Donato et al. 


Grandr & Falumbo 


2001 


12004 1 discuss the detec 


tion of a narrow FeKu line in the flat-spectrum radio quasar 3C273. 
In this case the thermal (unbeamed) and the beamed jet emission could 
be disentangled due to source variability, suggesting the presence of an 
underlying Seyfert-like component jSoldi et S? 2008 1 . 
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lOstorero et al.||7010[ |Kino et ar]|2013[ [Migliori et al.||2014| l 
and first source candidates have been discussed, e.g., 4C +55.17 
( McConville et al.||2011 1 with similar broadband properties as 
PMN J1603-4904. If confirmed as a y-ray loud young radio 
galaxy, then PMN J1603-4904 would be a well-suited object to 
investigate the origin of high-energy emission. 


5. Conclusions 


We have presented new X-ray observations of 
PMN J1603-4904, which further challenge its previous 
classification as a blazar and strongly suggest that this jet system 
is seen at a larger inclination angle to the line of sight. Its 
X-ray spectrum can be modeled by an absorbed power law and 
a Gaussian emission line at 5.44 + 0.05 keV, and can be best 
explained as emission from a non-blazar source. These X-ray 
properties are consistent with a young or evolved radio galaxy. 
We interpret the significantly detected spectral line as a neutral 
Fe Ka line with rest frame energy 6.4 keV. This identification 
results in the first redshift measurement for this source of 
z = 0.18 ± 0.01, i.e., a luminosity distance of 872 + 54Mpc. 

PM N J1603-4904 is a ssociated with a hard spectrum y-ray 
source (Nolan et al. 20121 1^ If it is indeed seen at a larger incli¬ 
nation angle where the observed emission is less intensified by 
Doppler boosting than in blazars, it adds to the class of so-called 
misaligned sources rarely detected at y-ray energies ( |Abdo et ak] 
2010^, with a remarkably high y-ray luminosity. 


Further multi wavelength observations are required to con¬ 
firm PMN J1603-4904 as a y-ray loud young radio galaxy: 
VLBI monitoring will allow us to check for proper motion in op¬ 
posite directions with respect to the core. Low radio frequency 
observations could help to determine a putative peak at < 1 GHz, 
thus increasing confidence in its identification as a GPS source, 
and to constrain the extended radio emission. 
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Table 2. Typical high-energy spectral parameters for different source classes compared to results for PMN J1603-4904. 
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